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Abstract In addition to the ordinary bulk higher equations of motion in the boundary 
^ ! version of the Liouville conformal field theory, an infinite set of relations containing the 
Q\ | boundary operators is found. These equations are in one-to-one correspondence with the 
| singular representations of the Virasoro algebra. We comment on the possible applications 
• ■ in the context of minimal boundary Liouville gravity. 

G\ '. 

O 1 Introduction 

^ \ Alexey Zamolodchikov showed in [1] that a special set of relations holds for quantum op- 
5_i ■ erators in the Liouville conformal field theory (LFT). They are parameterized by pairs of 
positive integers (m, n) related to the degenerate representation of the Virasoro algebra. In 
the classical limit, they represent the "higher equations of motion" (HEMs), because the 
first one (1, 1) coincides with the usual Liouville equation of motion. These equations relate 
different basic LFT primary fields V a (x). The equations are derived based on the two main 
postulates of LFT. The first is that all singular vectors vanish in the representations built 
on the exponential fields; this is a quantum version of the relations in the classical LFT. The 
second basically states that the set of exponential fields {V a } (with complex a allowed) in 
some sense covers the whole variety of primary fields in LFT. 

The higher equations turned to be useful in the context of minimal Liouville gravity 
(see [2] for the terminology). In particular, they were used in [3-5] to derive the general 
four-point correlation functions with one degenerate matter field. It is very likely that HEMs 
are potentially important in the general program of explicitly constructing the complete set 
of correlation functions in the minimal Liouville gravity. 
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Our purpose here is to reveal an additional set of relations for quantum operators in LFT 
with a conformal boundary (BLFT). In the next section, we collect the main facts about 
the BLFT. In section 3, we derive the boundary version of HEMs. Possible applications 
are considered in section 4. The definitions of the special functions and also some explicit 
calculations omitted in the main text are presented in the appendices. 



2 Boundary Liouville field theory 

We consider the Liouville conformal field theory on a domain T with a boundary dT. In a 
general background, the action of the BLFT is [6] 

Aound = ^J v [9 ab d a ^ + QR(j) + 47r/ie 2 ^] y/g#x + J (^-<j> + /i B e^ y/gidx . (2.1) 

ar 

The first term describes the theory in the bulk. Here, R is the scalar curvature associated 
with the background metric g, and fi is the bulk cosmological constant. The background 
charge Q = b + 1/b determines the central charge of the theory 

c L = l + 6Q 2 . (2.2) 

A conformally invariant boundary condition [7] is introduced through the boundary inter- 
action [6] given by the second term in (12. ip . Here, g\ is the induced boundary value of the 
metric, K is the geodesic curvature of the boundary, and Hb is the boundary cosmological 
constant. We always imply the upper half-plane geometry in what follows. In the bulk, the 
holomorphic component of the stress tensor has the form 

T(z) = -(<90) 2 + Q<9 2 0, (2.3) 

where z = x + iy. Because of the boundary condition on the Liouville field, the boundary 
value of the classical stress tensor is 

Td(x) = ~(pl + \p xx + 7r& 2 (7r/4& 2 - fie*, (2.4) 

where (p = 2b<p. This is equivalent to the classical equation 

+ Tcl ) e ~" /4 = ~ ( 2 - 5 ) 

relating two boundary exponents. This equation is an example of the classical limit of the 
quantum relations between the boundary primary operators, which we derive below. 

We let V a (z,z) denote the bulk primary fields. These fields have the conformal weight 
A a = a(Q — a). The structure constant C(a 3 , a 2 , a%) related to the three bulk primaries for 
generic values of a i: i — 1,2, 3, is given [8,9] by the expression 



-CHI —Otn — Cc-J 



C(a 3 ,a 2 ,ai) = n^(b 2 )b 2 2b2 

T T b (2a 1 )T b (2a 2 )T,(2a 3 ) 



T b («i + a 2 + a 3 - Q)Y b (a 1 + a 2 - a 3 )T b (o;i + a 3 - a 2 )T b (a 2 + a 3 - oti) ' 



(2.6) 

W3) J- -r « 3 — tx 2 ) i b ytx 2 -r u!3 — «lj 
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where the special function T b and also the other special functions used in the main text are 
defined in Appendix A. 

The boundary operators B^ 2CT1 (x) have the conformal weight = f3{Q — (3) and are 
labeled by the two indices o\ and <r 2 related to the left and the right cosmological constants 
fi Bl and /j,b 2 by 

cos (2irb(a - = ^^sm(nb 2 ). (2.7) 

The observables depend on the scale-invariant ratios fJ>i/ii B - For example, the correlation 
function of n bulk operators V ai . . . V an and of m boundary operators B^ 2 . . . B a ^ ai scales 



as 



\ fi H V J 

where F is some scaling function. 

To characterize the LFT on the upper half-plane, we must know [6] some other structure 
constants in addition to the bulk three-point function C(a±, a: 2 , 0:3): 

1. the bulk one-point function [6] 

\z — z\ 

2. the boundary two-point function [6] 

S(p 1 ,a 2 ,a 1 )5(p 2 -p 1 ) 



(b;^(x)b;p(o))- 

3. the bulk-boundary two-point function [10] 

<v„(*,*)fly(*)>- R[aJ3M 



and 

4. the boundary three-point function [11] 

£l(CT3<T2Cri)/?3 

/ D 0-10-3 / \Dff3ff2/ T \R0 2 oi/ \\ _ ftffl 

\ D Q-P3\ X Z) D I32 K x V n Pi y x D/ - 1 1A1+A2— A3 I 1A2+A3— Ai I l A 3 +A 1 -A 2 - 

\ x 2l\ \ x 32\ p3l| 

We here give the boundary structure constant explicitly because we use it in the following 
sections: 

c (030 2 0i)/3 3 = ( w ( 6 2) 6 2-26 2 )a(A-A-ft) 

r 6 (2g - a - & - /3 3 )r 6 (/3 2 + & - A)r 6 (Q + fc-Pi- fo)T b (Q + p 3 - Pi - ft) 



X- 



X 



r b (2ft - Q)r 6 (Q - 2/? 2 )r 6 (g - 2p 1 )r b (Q) 

Sb(Jh + 0-1 - o-3)Sb(Q + P3 - 0-3 - 0-1) 

S&(/3 2 + 0"2 - °3)S b {Q + P2 - - 0"2) 

ioo 

X « i S S^V, + s)S 6 (V 2 + s)S b (V 3 + s)S 6 (Q + s) ' 1 ' 

—ioo 
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The coefficients Ui and Vi, i = 1, . . . , 4 are 



[/ x = 01 + <j 2 - Vi 

U 2 = Q -a 1 + a 2 - Pi, V 2 

U 3 = + cr 2 - cr 3 , V3 
t/4 = Q - + 0"2 - 0"3- 



Q + ct 2 - o"3 - A + /3 3 , 
2Q + (j 2 -a3-/3i-/53, 
2a 2 , 



3 Higher equations of motion 

Before discussing the boundary case, we briefly recall the reasoning leading to HEMs in the 
LFT without boundary [1] . The degenerate primary field V m ^ n appears for the Kac values [12] 
of the conformal dimension A m n related to the parameter 

Q (mb- 1 + nb) 

«m,n = — ~ • (3.1) 

We let D m ^ n denote the singular vector creating operators. It was proved in [1] that the 
primary field 

D m ,nD rn ^ n V mn (3-2) 
can be constructed, where the degenerate logarithmic field V' a is defined as 



y, 1 9V a 



2 da 



(3.3) 



Field (13.21) has the conformal dimension A m , n + mn. It should be identified with the primary 
field V m - n of the same dimension. The operator-valued relation 

Dm,nD mjr iV m n B m n V m — n (3.4) 

then holds. The coefficients B m>n are defined explicitly: 

m— 1 
n-1 

B mtn =(7c^(b 2 )) n b 1+2n - 2m ^m-nb 2 ) H (Ib^ + kb). (3.5) 

k=l— n 
1=1— m 

(fc,o^(o,o) 

To formulate the boundary analogue of the bulk HEMs, i.e., the relation that should 
relate the boundary operators, we first note that the method described above is no longer 
applicable. Indeed, for the bulk construction we need both left and right singular vector 
creating operators, while we have only one Virasoro algebra on the boundary. Instead we 
consider the action of D m ^ n on the primary boundary field B^ s ^ with arbitrary values of the 
boundary parameters. The analysis [6] of classical limit (12. 5p shows that this field should 
not vanish in the general case. On the other hand, it follows for purely algebraic reasons 
that the field D mjn B^ has the properties of the primary field regardless of the value of the 
boundary cosmological parameters. Taking the main assumption of LFT into account, i.e., 
there exists only one primary field of a given conformal dimension, we must identify 

/-) f>S lS2 _ TfS 1 S 2 f>S 1 S 2 (O C\ 
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up to a numerical constant, where we have the primary boundary field of the same conformal 
dimension in the RHS. This operator- valued relation assumes the corresponding relations 
between correlation functions if one of the fields is the subject of (I3.6p . In principle, the 
consistency of this statement must be verified for arbitrary correlation functions. As usual, 
it suffices to consider the consequence of (13.61) for the three-point functions or, equivalently, 
for the structure constants of the boundary operator product expansion 

(B^(0)D m>n B^(x)B s ^(oo)) = (B^(0)B s ^)B^(oo)) . (3.7) 

We use this relation to define the coefficients 

m > n ~ (B^(0)B^ n (x)B^(oc)) • 1 ' ] 

We note that a rather nontrivial consequence of (13.61) is that the ratio of two correlation 
functions (13. 8ft depends neither on the conformal dimension of the other two fields (3\ and 
/?2 nor on the cosmological parameter S3 of a boundary segment not directly connected to 
Bm S n( x )- Before considering the general situation, we test this idea in the case where the 
screening calculations allow avoiding the complicated special functions in expression (12. 8p 
for the general three-point boundary correlation function. 



3.1 Screening calculations 

In this section, we apply (13.61) in the case where the three-point correlation functions can be 
computed perturbatively as Coulomb gas integrals [13,14]. We recall that if the conformal 
parameters ccj and /3j of the correlation function (Vi • • • B\ • • • ) satisfy the screening relation 
E a i + EA = Q ~~ n b> then this correlation function has a pole [15], and the residue is 
calculated using the perturbation theory in /i and hb- Because of the total charge balance 
condition, only a finite number of terms in the series have nonzero values. We consider 

D lt2 B^ 2 = K{]° 2 B{] s _\. (3.9) 

Taking into account that ai j2 = —6/2 and 01,-2 — 36/2, we chose the other two fields such 
that the screening relation 

(5--(0)^i, 2 n?2(^)^-3 b /2( 00 )) = K tT (^ 3Sl (0)^ / 1(^)^_ 3f)/2 (oo)) (3.10) 

is satisfied. The action of D 12 reduces to a factor that is known explicitly [1]. The total 
charge balance is performed for the correlation function in the RHS, and hence 

K S ,T = 2(1 - 26/3)(l - 26/3 - 6 2 ) (B?«{0)B™{l)B™_ap{oo)) " ( 3 - n ) 

A nontrivial check of the general statement should be that the dependence on (3 and /13 in 
the RHS of (13. lip cancels. The volume screening contribution is related to the interaction 
— \i J d 2 ze 2b< ^, while the boundary contribution comes in the second order and requires two 
boundary screenings e^, 

(s--(o)5i i b 7 2 (i)^-3 & / 2 (°°)) = -a* / ^^(^ ) ^ 3S H0)^; 2 (i)^-3 b/2 ( oo )) 

lmz>0 

+ E^f L / / ^i^e^e^)^ (3.12) 
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where the contours are defined as C\ = (— oo,0), C 2 = (0,1), and C 3 = (l,oo) and ^ are 
the corresponding values of the boundary cosmological constant. The explicit expressions 
for the free theory correlation functions are 

(e 2 ^)S--( )S-7 2 (l)^-_ 36/2(oo)>0 = \ Z \-W |! _ z f \ z _ z{ -^ (3 . 13) 

and 

(e 60(xi) e 6 ^-)^-(O)5-7 2 (l)^_ 36/2(oo))o = 

I I —26/3 I I —26/3 I -i i6 2 I-, 16 2 i i-26 2 /, 1 

— Fl| F2| 1 1 — rCl ] |1 — X2\ \X\—X2\ ■ (3.14) 

We introduce the notation 

I(A,B,C)= J d 2 z\zz\ A \(l -z)(l -z)\ B \z-z\ c , 

lmz>0 

hj(A : B, C) = J c J c dx\dx2 \x\\ A \x2\ A |1 — x\\ B |1 — x 2 |' B \x\ — X2\ C ■ (3.15) 
We can write 

(B£«$)B?$ 2 tx)B%« _ 36/2 (oo)> = 

= -tf(-2bp, b\ -2b 2 ) + J2 ^Iv(-2bp, b\ -2b 2 ). (3.16) 

All integrations can be performed explicitly. Using the results in Appendices C and D, we 
obtain 

/ = — \r sin(— ) sin(TrA) sm(nB) sin(7r(A + B + C))J(A, B, C), 
2tt 6 2 

1 llC C 

In = 7 sin — smir(A + B + — ) sin ir(A + B + C)J(A,B,C), 

it A 2 2 

1 C C 

In = sinirC simr{A + B + -)sin7r(A+ —)J(A,B,C), 

I13 = sin nC sin n(A + B + ^)smn(B + J(A,B,C), (3.17) 

2ii' i 2 2 

\ ixC C 

I22 = o sin -— sin7r(A + — ) sirtixAJ(A, B, C), 

it A 2 2 

1 C C 

hs = - ^3 s in ttC sin tt(A + -) sin tt{B + -) J(A, B, C), 

1 llC C 

I33 = » sin — sin7r(B + — ) sin %BJ(A,B,C), 

7T° Z Z 



where 



j(A c) = r(A + i)t(b + i)r(c + i)r(-c/2)r(s + c/2 + 1) x 

T(-A -B-C- 1)T(-A -B-C/2- 1)T{A + C/2 + 1). (3.18) 
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Summing these contributions, we obtain 



(^(Ojfl^B^Joo)) 



r sin — 

2vr 3 2 



&/2V^-"Q-/3-3&/2^ 
7I"C 

- fi sin(7rA) sin(7T-B) sin(7r(v4 + B + C))+ 



+ ^ sin(7r(A + B + C/2)) sin(7r(A + 5 + C))+ 

+ ^ sin(7r(A + C/2)) sin(vrA) + /i 2 sin(7r(S + C/2)) sin^S)- (3.19) 

71~C 

- 2/ii/i 2 cos — sin(7r(A + B + C/2)) sin(7r(A + C/2))- 

7rC 

- 2/1x112 cos — sin(7r(A + B + C/2)) sm{ix{B + C/2)) + 



7rC 

+ 2/i 2/ u 3 cos — sin(vr(A + C/2)) sin(7r(5 + C/2)) 



J(A,B,C), 



where A = —2b(3, B = b 2 , and C = —2b 2 . It is sufficiently remarkable that because of the 
relation C = —25, the result is independent of /z 3 , 



i 



/3-3b/2 v 



27T 3 



/i sin(7r_B) + + n\ — 2/i!/x 2 cos(7r5) 



x sin(TrA) sia{irB) sin(7r(A - 5)) J(A, B, C). (3.20) 



With the FZZ parameterization 



2 cosh 7T0Sj 

Pi = P — : W 

sin irb z 

correlation function ( 13.121) is 

(5; 3S1 (0)5-(,)B-_ 3W2 (»)) 



cosh nbsj cosh Trbs,-, 



sin 7r6 2 



(3.21) 



b/ 2 VV-"Q-/3-3&/2^ 

1 i? + zsi — zs 2 -B — «Si + zs 2 -B + zsi + iS2 B — isi — zs 2 
sin sin sin sin 



2tt 3 ~~" 2 2 2 

x sin(TrA) sin(TLB) sin(7r(A - B))J(A, B, C). 



(3.22) 



It can be easily verified that /3-dependence of 2 also vanishes. Finally, 

b - i(s 2 - sx) 



S\S 2 



71 



. ,b + i(s 1 + s 2 ) . ,6-i(si + s 2 ) . ,& + H s 2-si) . 

sin 7ro sin 7to sin 7ro sin ixb- 

2 2 

A similar calculation in the case (1,1) gives 

/4/i 7 (o 2 ^ 1/2 



(3.23) 



^1,1 



1*2 



V 7r 



1(1 — sin irb sin irb 



(3.24) 



In the next section, we show that (13. 23j) and (I3.24[) are generalized for the general case (m, n). 
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3.2 General three-point analysis 

Operator-valued relation (13.61) means that the equality ( 13 .TP holds for the general three- 
point correlation functions. In terms of the boundary structure constants (12.81) this gives the 
following expression for the coefficients 



^m,n _ Pm,n(Q ~ @3 ~ Pl)Pm,n\Pz ~ Pi 



(s3«2Sl)/33 



&m.,n$l 



(s 3 S 2 Sl)/33 ' 



a 



(3.25) 



m,,n-\-nb 



.01 



where P m ,n is the fusion polynomial 

Pm. .n ( X I 



(3.26) 



fc=l— n:2:n— 1 
l=l—m:2:m—l 



and A; 5 fc = (Zfe 1 + kb)/2. To calculate the ratio, we use the following generalizations of the 
shift relations for the S& and r& functions presented in Appendix A: 



n-l 



S&(a; + nb) = 2 n JJ sin 7rfe(a; + fen) • 5* ( 



(3.27) 



fc=0 



and 



r b (x + nb) = { Z!J^„, - — r b (xV 



r + **)] 



(3.28) 



It is convenient to split the ratio into three parts. The first contains the ratio of the integral 
parts in expression (12. 8p . Here, the diference comes from two Sb functions in the integral. 
Keeping in mind that Sb(Q — x) = 1/S b (x), we obtain 

Sb(Pmn + nb + (J 2 - 0"3 + s)S b (Q - P mn ~ Tib + <T 2 - <J 3 + s) = 



sm 



n-l 

n 

fc=o sm 



nb 



(l+2fc-n)b 



+ (7 2 - (7 3 + S + 



1— m , 



-7T 



— nb 



(l+2fc-n)b 



1— m 



7T 



2 + CT 2 - CT 3 + S 1 + 2 
n + CT 2 - CT 3 + s)£&(Q - Pmn + <T 2 - CT 3 + s) 

(-l) m "^(/5 mn + a 2 - a 3 + S )S fe (Q - (3 mn +a 2 - a 3 + s). 



(3.29) 



Hence, the ratio of the integrals gives just (— l) mn . The second part of the ratio is the part 
containing the S b functions in the prefactor in front of integral ( 12. 81) . It contributes 

S b (Pmn + nb + a 2 - a 3 )S b (Q + (3 mn + nb - a 2 - cr 3 ) 



S b (f3 mn + cr 2 - a 3 )S b (Q + P mn - c 2 

n-l 

2 2n JJsm 7Tl> \ ; a 2 - rr ?l 

k=0 



^3) 



sm 



nb 
Tib 



1 



m 



(3 + 2k-n)b 



a 2 - a 3 



2 

3 — m 



-IT 



-7T 



(3.30) 
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to the ratio. In the FZZ parameterization (I3.2ip parameter a = Q/2 + is/2 and we have 



r ' (1 - m)^ 1 + (1 + 2k - n)b + i(s 2 - s s 
no 



2 2n J] si 



sin 

fc=0 

• sin 



Tib- 



2 

1 - m)& _1 + (1 + 2k - n)b + i(s 2 + s 3 ) 



(3.31) 



2 

The third part of the ratio, the part containing the Tj, functions in the prefactor, contributes 

M _ T b( 2 Q -Pi- Ann ~ Ps^bjPmn + fa - gQ 

T b (2Q -f3 x - (3 mn -nb- f3 3 )T b ((3 mn + nb + (3 3 - fa) 
T b (Q + (3 mn -Pi- (3 3 )V b (Q + (3 3 - (3mn - f3i)T b (Q - 2{3 mn - 2nb) 
' T b (Q + (3 mn + nb - ft - (3 3 )T b (Q + (3 3 - (3 mn - nb - (3x)V b (Q - 2(3 mn ) ' 

The dependence of M on f3\ and f3 3 comes from two factors of the form 

T b ((3 mn + (3)T b (Q - (3 mn + (3) = b- nb ^+ nb -Q)b~ mn 
FbiPmn + nb- p)T b (Q - (3 mn - nb + (3) Pmn(P) 

with (3 equal to either (3 3 — (3\ or Q — (3$ — fti, and hence 



(3.32) 



(3.33) 



M = f)- 2nb ^l3 m n+nb-Q) -2mn \ ^b(Q ZPmn 2n6) ^ 

Pmn(^3 ~ ^l)Pmn(^3 + -^l) ^b(Q ~ ^Pmn) 



Combining all together, we obtain 

2n-l 



= {-l) mn (^Ll^l\ 2 fe 2„(l-m) "Q r(m _ ( n _ (3>35) 
V ^ / i,_n 



fc=0 

n-1 



n. , (1 - m)b~ l + (1 + 2k - n)b + i(s\ + s 2 ) . (l-m)b- 1 + (l + 2k-n)b + i(s 2 -s 1 ) 
Sin 7T0 SUl 7T0 - 
9 9 



k=0 



It is easy to see that the degenerate field B^'* 2 has a vanishing singular vector and the 
truncated operator product expansion if si ± s 2 = 2i\k^ r with k — 1 — n, 3 — n, . . . ,n — 1 and 
r = 1 — m, 3 — m, . . . , m — 1. This result generalize (1, 2) fusion rules suggested in [6]. 

4 Application in minimal gravity 

One possible application of the BHEM in a physical context is for constructing the corre- 
lation functions of physical fields in boundary minimal Liouville gravity (BMLG). This is 
an alternative description of non-critical open string [16] propagating in "low-dimensional" 
space-time. In BMLG, the gravity is induced by one of the minimal CFT models, and it is 
expected to be exactly solvable. This is confirmed by the fact that the alternative approach 
of matrix models to 2D gravity (see, e.g., [17] for a review) provides explicit expressions for 
many observables. The comparison of some open string amplitudes, derived both in matrix 
model framework and using worldsheet description, was performed recently in ref. [18]. The 
BMLG theory consists of the matter, Liouville, and ghost sectors, which do not interact 
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except through the conformal anomaly and through the constriction of the physical fields. 
The matter central charge is defined by the central charge balance condition 



c M + c L = 26. (4.1) 

The physical fields are defined in the framework of BRST quantization as cohomologies with 
respect to the BRST charge 

Q = j (C(T L + T M ) + CdCB) (4.2) 

Here, B and C are ghost fields of the respective spins 2 and -l(we use the uppercase letters 
for the ghosts here in order not to confuse B with the Liouville parameter b). 

It is the specific property of MLG that in the construction of the physical fields of 
nonzero ghost number all matter fields are "dressed" by Liouville exponentials in the LHS 
of the HEMs. For example, there exist ghost number-1 basic boundary physical fields of the 
form 

and 

* m,n ^m,~n- V / 

Here, the parameters ai,a 2 and s±,s 2 correspond to the conformal boundary conditions 
to the left/to the right from the operator insertion in the respective matter and Liouville 
sectors. Because of the anomaly of the ghost current, the MLG correlation function of any 
number iV of fields must be of the form [19] 

N 

G N = J] / dx i (W 1 (x 1 )W 2 (x 2 )W 3 (x 3 )U 4 (x 4 ) . . . U n (x n ))mG: (4.5) 

i=4 •* 



where (. . -) MG denotes the joint correlation function of matter, Liouville, and ghosts. 

Another important class of boundary physical fields (ghost number 0) is the boundary 
ground ring. The general form of the elements of the boundary ground ring [20,21] (also 
see [22] for more recent developments) is 

O m ,n H mn ty m,nVm,n- (4-6) 

Here, H m ^ n are operators of level mn and ghost number constructed from the Virasoro 
generators L^t\ and ghosts. It can be shown analogously to [23] that the boundary 

higher equations lead to the following important relation between the two types of physical 
fields introduced above: 

QQSl,s 2 = Ksuszyysusz U 7) 

-<-' v ^m,n m,n "ran ' V J 

In particular, this means that to construct the boundary ground ring element O^^ 2 , the 
fusion relations for the cosmological constants s 1 and s 2 in the Liouville sector should be 
satisfied, K^£ 2 — 0- Another consequence of this relation is that if the fusion rules for s 1 
and s 2 are not satisfied, then the field W^£ 2 seems exact, and the correlation functions of 
this field are naively equal to zero. This statement should be checked more carefully. 
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Taking the commutation relations {B_i, Q} = d into account, we can straightforwardly 
verify that 

ufar^ = (d - cb_o o£r |si,S2) (4-8) 

This relation allows performing every (integrated) insertion of the particular operator Uk(xk) 
explicitly, integrating by parts. The integral (14. 5 ft can thus be reduced to boundary terms, 
which are in principle defined by the operator product expansions of the ground ring elements 
and the basic boundary fields W a . 
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Appendix A. Special functions 

Here we collect the definitions and some properties of the special functions r 6 , T b and S b . 
The Double Gamma function introduced by Barnes [24] is defined as 



f, oo 

a 



logr 2 (s|u; 1 ,^ 2 )= — ^ (s + + n 2 u 2 ) 



dt 

The Barnes Gamma function is defined as 



ni,n2=0 / £_ 



The function V b (x) satisfies the functional relations 



T b (x + b) 



T b {x + 1/6) 



2nb bx - 



T(bx) 



-Th(x) 



T(x/b) 



r b (x) 



and is a meromorphic function of x, whose poles are located at x — —nb — mb 1 ,n,m G N. 
The integral representation converges for < Rex: 



logr 6 (x) 



00 dt 
T 



—xt 

e — e 



-Qt/2 



,-bt 



)(l-e-'/ b ) 



(Q/2-x) 2 c _ f _ Q/2-x 
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The Barnes Sine function is defined as 



S b (x) = Tbix) 



T b (Q-x) 



(A.2) 



The shift relations are 



S b (x + b) = 2sm(nbx)S b (x), 
S b (x + 1/6) = 2sm(nx/b)S b (x), 

and S b (x) is a meromorphic function of x, whose poles are located at x — —nb — mb~ x ,n, m e 
N, and whose zeros are located at x = Q + nb + mb~ l ,n, m e N. The integral representation 
converges in the strip < Rex < Q: 



f 00 dt 

\ogS b {x) = j 



sinh(§ - x)t (Q - 2x) 



2sinh(f)sinh(^) 



Finally, the upsilon function is defined as 

T h {x)- x = V h {x)V h {Q-x). 

The functional relations are 

r b (x + b) = J} bx l ^ -^(x), 



(A.3) 



r b {x + 1/6) 



r(i - bx) 
r(i -x/b) 



b 2x ' b - l T b {x), 



and Tfe(x) is an entire function of x whose zeros are located at x — —nb — mb 1 and 
x = Q+nb+mb' 1 , n, m E N. The integral representation converges in the strip < Rex < Q: 



logT, 



Jo 



dt 



Q 



sinh 2 (§-x)| 
/ sinhfsinh^ 



Appendix B. Contour calculus. Bulk term 

We consider the integral 

1=1 d 2 z{zz) A [(1 - z)(l - z)] B \z - z\ c . 

Jlmz>0 

It is convenient to extend the integration area to be the whole complex plane: 

r\C roc 

^ I , / o o\ A r / - \o 21 B C 



(B.l) 



1 + e 



ivrC 



/oo 
dxdy(x 2 + y 2 ) A [(l- x y + y*] B y 
■oo 



(B.2) 



For practical computations, it is convenient to decompose this integral into a sum of holo- 
morphic and antiholomorphic parts using the way proposed in [25,26]. Performing the Wick 
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rotation y — > iy(l — lie) (with infinitely small e > 0) and introducing the new variables 
u = x — y and v = x + y, we can write the integral in the factored form 



/oo />oo 
dv / <i-u(t> — — u)t) A {u + i(v — u)e) J 
-OO J —CO 



2(1 + e- c )^..-.. 

(1 - u - i(u - w)e) B (l - u + i(u - w)e) B (w - u - i(u - u)e) c . (B.3) 



There are three branch points u = —iev, u — 1 — — l)e, and u = v + ie in the w plane. By 
deforming the integration contour, we can easily demonstrate that nontrivial contributions 
come only from two domains of integration, namely, when < v < 1 and v > 1, 



2(1 + e^ c ) 

where 



[h + h], (B.4) 



duu A {\ -u) B {u-v) c 



h = J* dvv A {l - v f^e i7rC jf * duu A (l - u) B (u - v f + e- mC 

/»oo n 1 

+ e tn(C + B) / duu A( u _ X )B( U _ + e -«r(C+B) / duil A (u — l)^(w — f)'"' 

</ 1 J 00 

/oo f />oo 

+ e l7r( - B - c) y"^ A ( M -l) B ( M -t;) c |. (B.5) 



In a more compact form, 



J 1 = 2*sin7rC [ dvv 2A+c+1 (l-v) B f duu A (l - u) c (I - vu) B + 
Jo Jo 

+ 2isimi(C + B) f dvv A {l-v) B f duu^ A - B ~ c '\l - u) B (l - vuf , 
Jo Jo 



h = 2i sin ttC f dvv- 2A ~ 2B ~ c -\l - v) B f duu A ' B - c - 2 {\ - uf{\ - vu) B . (B.6) 
Jo Jo 



We obtain the expression 

I = — cos — 
2 V 2 



sin vr(C + B)I (A, B, -A - B - C - 2, B, C) + 
sin nCI (2A + C + 1,B,A, C, B) + 
sin nCI (-2A - 2B - C - 3, B, -A - B - C - 2, C, B) 



(B.7) 



where 

io(a,/3,7,£,e) = I dvv a {l-vf! duu\l - u) 5 (l - vu)\ (B.8) 
Jo Jo 
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We now use standard formulas for the integral representations of the higher hypergeometric 
functions 

I {a,/3, 7, S, e) = B{ 1 + 1,5+1) / dvv a (l - vf 2 F 1 (-e, 7 + 1, 7 + 6 + 2; x) = 

Jo 

5(7 + 1,5 + l)B(a + 1,(3 + l) 3 F 2 (-e, 7 + 1, a + 1; 7 + 8 + 2, a + /3 + 2; 1), (B.9) 

where B(x,y) is the beta function. In our cases, the integral J appears at special values of 
the parameters and can be expressed in terms of gamma functions because of the formula 

r(l + q/2)r(l + a - 6)r(l + a - c)r(l + a/2 - b - c) 

3 F 2 (a, 6, c; a - b + 1, a - c + 1 1) = — — - — — — -. 

T(l + a)r(l + a/2 — 0)1(1 + a/2 — c)T(l + a — b — c) 

(B.10) 

Thus, 

I (A,B,-A-B -C-2,B,C) = sin7rCsin7r(A + B + — ) smn(A + — ) • J(A,S,C7), 

27r J 2 2 

I (2A + C + 1, B, A, C, B) = sin sin n(A + B + -) smn(A + B + C) ■ J(A, B, C), 

2ix 6 2 2 

(B.H) 

1 ttC 1 C 1 

h(-2A - 2B - C - 3, B, -A - B - C - 2, C, B) = sin sin tt(A + -) sin nA ■ J(A, B, 

2t: 6 2 2 

where 

J(A, B, C) = Y(A + 1)Y(B + 1)Y(C + l)V{--)V(-A -B-C- 1)T(-A -B---1) 

xT(B + ^ + l)T(A + ^ + l) 

(B.12) 



Substituting (lETTTD in (TBTTj) gives 



1 / ttC x 



.cos J sin nA sin irB sin nC sin ir(A + B + C)J(A,B,C). (B.13) 

Appendix C. Contour calculus. Boundary terms 

Here, we give the details of the calculations of the integrals that appear when calculating the 
coefficient K12. The parameterization of the boundary is chosen as follows. The intervals 
(—oo,0), (0, 1), and (l,oo) are respectively denoted by C\, C 2 , and C3. The corresponding 
cosmological parameters are S3, si, and s 2 - We are interested in the integrals 

hj = / dxi / c?X2|a;i| A |a;2| A |l — — X2| B |xi — X2| c '. (C.l) 

Ja Jd 
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With = Iji taken into account, it suffices to consider i < j. Below, we often use the result 
for the two-dimensional Selberg integral 

dx! [ dx2\x 1 \v- 1 \x 2 \^ 1 \l -aril" -1 !! - x 2 \ v ~ 1 \x 1 - x 2 \ 2g = 



T(g)r(fi + u + g)T( f i + u + 2gy 1 ' } 

The first integral has the form of the Selberg integral and can be written as 

1 ttC C 
in = ? sin — sin n(A + B + — ) sin n(A + B + C)N. (C.3) 

The integral ii2 reduces to the function Iq introduced in (IB. 91) (in the previous appendix) 
after x 2 — ► 1/^2, 

i 12 = i (A, B,-A-B-C-2,B,C). (C.4) 

To calculate ii3, we consider the contour integral in the complex plane x 2 such that the 
contour goes along the real axis above the branch points 0, x\, and 1. The contour can be 
moved to infinity, which gives the relation 

pi px\ pi pi pi poo 

J 13 + e -™A / / +e ~in(A + C) / / +e -i„(A+C+B) =0 (C.5) 

JO J0 Jo Jxi Jo Jl 



or 

"1 px\ 

Ii3 + e~™ A 

io Jo 





pi px\ - 


+e ^in(A+C) 


hi- / 




Jo Jo . 



+ e -MA+C+B)j i2 = Q (C6) 



In the interval (0,xi), the integral coincides with Iq{2A + C + 1, B,A,C, B). Taking into 
account that all integrals are real, we obtain 

ii 3 = - cos n(A + (7) in - cos ix(A + B + C)ii2+ 
+ [cosvr(A + C) - cos 7vA}I (2A + C + 1,B,A,C,B). (C.7) 

Integral I 22 reduces to the Selberg integral after x 2 — ► l/x 2 and x 2 — ► l/x 2 , 

1 7rC C 

J 22 = -sin sin7r(A+ — ) sinvrA J(A, 5, C). (C.8) 

ix 6 2 2 

To calculate i 2 3, we do the same trick as for i 13 . The relations 

POO pi POO PX\ POO POO 

i 23 + e~ inA / / +e-™ {A+B) / / +e -^(A+c + B) / / = o (C.9) 

Jl Jo Jl Jl Jl Jx 1 

and 



OO POO POO POO 



I 23 + e-^ A I 12 + e-^ A+B ^[I 22 - / ] + e -**(.A+c+B) / =0 (ai0) 

Jl Jx 1 Jl Jx 1 
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hold. In the interval (x\, oo), the integral coincides with I (—2A — 2B — C — 3, B, —A — B — 
C-2,C,B). Hence, 



I 23 = - cos irAI 12 - cos tt(A + B)I 22 + [cos 7f(A + B) - cos tt(A + B + C)\ 

■h{-2A - 2B - C - 3, B, -A - B - C - 2, C, B). 



(C.11) 



The integral -Z33 is found from the relation 



r xj 



+e 



-inC 



— oo J — oo 



rO 



+e 



-iTT(C+A) 



OO J Xl 



/-I 



oo </0 



+e 



-i7r(A+C+B) 



roo 



oo J 1 



0. 



Taking into account that 



we obtain 



PX 1 



—00 j —00 



rO 



' 33 



00 »/ Xl 



1 + e 



-iTrC 



j 33 + e -^ c+A h l3 + e -^(^+c+s) J: 



23 



(C.12) 



(C.13) 



(C.14) 



and 



'38 



COS 



c c 

COS7T(A + -)/ 13 + COS7T(A + 5 + -)/ 23 



(C.15) 



In the calculation of the boundary integrals, we have thus obtained the expressions 



In 
I12 

/l3 
I22 
^23 
I33 



1 7rC C 

-—sin sin n(A + B + —) sinvr(A + B + C) J(A, 5, C), 

ix 6 2 2 

1 C C 

-— sin7rC sin7r(,4 + B + — ) shi7r(A + — ) J(A,B,C), 

1 C C 

sin rtC sin ir(A + B + — ) sin 77(5 + — ) J(A, B, C), 

2tx a 2 2 



7T" 



sin sin ir(A H ) sin nA J(A, B : C), 



1 C C 

sin tiC sin it(A + — ) sin ^(5 + — ) J(A,B,C), 

2^^ 2 2 

1 7rC C 

sin sin tc(B H ) sin 7TB J(A, fi, C). 

n s 2 2 



(C.16) 
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